AD-A198  579 


i  - 


mC.ElLt.Uiia 


AFWAL-TR-88-2007 


A  NUMERICAL  SIMULATION  OF  THE  FLOW  FIELD 
AND  HEAT  TRANSFER  IN  A  RECTANGULAR  PASSAGE 
WITH  A  TURBULENCE  PROMOTER 


Bryan  R.  Becker,  Ph.D.,  P.E. 
University  of  Missouri 
Kansas  City,  MO 


Richard  B.  Rivir,  Ph.D. 
Components  Branch 
Turbine  Engine  Division 


s 


DTIC 

ELECTEI 

AUG  0  5  1988 G 


% 


June  1988 


Final  Report  for  Period  June  1987  -  September  19S7 


Approved  for  public  release;  distribution  unlimited. 


AERO  PROPULSION  LABORATORY 

AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 

AIR  FORCE  SYSTEMS  COMMAND 

WRIGHT-PATTERSON  AIR  FORCE  BASE,  ( HIO  45433-Cbt3 


NOTICE 


WHEN  GOVERNMENT  DRAWINGS,  SPECIFICATION'S,  OR  OTHER  DATA  ARE 
USED  FOR  ANY  PURPOSE  OTHER  THAN  IN  CONNECTION  WITH  A  DEFINITELY 
GOVERNMENT-RELATED  PROCUREMENT,  THE  UNITED  STATES  GOVERNMENT 
INCURS  NO  RESPONSIBILITY  OP.  ANY  OBLIGATION  WHATSOEVER,  THE  FACT 
THAT  THE  GOVERNMENT  MA-Y  HAVE  FORMULATED  OP  IN  ANY  WAY  SUPPLIED 
THE  SAID  DRAWINGS,  SPECIFICATIONS,  OP.  OTHER  DATA,  IS  NOT  TO  BE 
REGARDED  BY  IMPLICATION,  OP  OTHERWISE  IN  ANY  MANNER  CONSTRUED,  AS 
LICENSING  THE  HOLDER,  OR  ANY  OTHER  PEPSON  OR  CORPORATION;  OR  AS 
CONVEYING  ANY  RIGHTS  OR  PERMISSION  TO  MANUFACTURE,  USE,  OR  SELL 
ANY  PATENTED  INVENTION  THAT  MAY  IN  ANY  WAY  BE  RELATED  THERETO. 

THIS  REPORT  HAS  BEEN  REVIEWED  BY  THE  OFFICE  OF  PUELIC 
AFFAIRS  CaSD/CPA)  AND  IS  RELEASABLE  TO  THE  NATIONAL  TECHNICAL 
INFORMATION  SERVICE  CnTIS) .  AT  NTIS,  IT  WILL  BE  AVAILABLE  TO  THE 
GENERAL  PUBLIC,  INCLUDING  FOREIGN  NATIONS. 

THIS  TECHNICAL  REPORT  HAS  BEEN  REVIEWED  AND  IS  APPROVED  FOF 
PUBLICATION . 


Project  Engineer  Chief,  Components  Branch 

FDR  THE  COMMANDER 


R.  E.  HENDERSON 
Deputy  for  Technology 
Turbine  Engine  Division 


IF  YOUR  ADDRESS  HAS  CHANGED.  IF  YOU  WISH  TO  BE  REMOVED  FROM  OUR  MAILING 
LIST  OR  IF  THE  ADDRESSEE  IS  NO  LONGER  EMPLOYED  BY  YOUR  ORGANIZATION  PLEASE 
NOTIFY  AFWAL/POTC  .  WRIGHT -PATTERSON  AFB.  OH  45433-  6563  .  TO  HELP  MAINTAIN 
A  CURRENT  MAILING  LIST. 

COPIES  OF  THIS  REPORT  SHOULO  NOT  BE  RETURNED  UNLESS  RETURN  IS  REQUIRED  BY 
SECURITY  CONSIDERATIONS.  CONTRACTUAL  OBLIGATIONS.  OR  NOTICE  ON  A  SPECIFIC 
DOCUMENT. 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No  0704-0188 


la  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED _ 

2a  SECURITY  CLASSIFICATION  AUTHORS 

N/A _ 

2b  DECLASSIFICATION /DOWNGRADING  iCHEDUl  F 


1b  RESTRICTIVE  MARKINGS 

N/A _ 

3  DISTRIBUTION /AVAILABILITY  OF  REPORT 

a r pmvf-rl  for  public-,  distribution 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

A  AL  -TR-88 -  2  00  7 

6a  NAME  OF  PERFORMING  ORGANIZATION  1 6b 


5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


Components  Branch 

6c.  ADDRESS  {City,  State,  and  ZIP  Code) 


6b  OFFICE  SYMBOL  1  7a.  NAME  OF  MONITORING  ORGANIZATION 
(If  applicable)  \ 


AFWAI/POTC 


7b  ADDRESS  (City,  State,  and  ZIP  Code) 


AFWAL/POTC 

Wright -Patter  son  AF3  011  45433-6563 


AFOSR/NA 


10  SOURCE  OF  FUNDING  NUMBERS 

PROGRAM  I  PROJECT  |1 

ELEMENT  NO  NO  I 

61102F  2307 


PROJECT 

TASK 

WORK  UNIT 

NO 

NO 

ACCESSION  NO 

2307 

S3 

1^ 

8a.  NAME  OF  FUNDING /SPONSORING  8b  OFFICE  SYMBOL  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

.  ORGANIZATION  (If  applicable) 

AFOSR  AFOSR/NA 

8c  ADDRESS  (City,  State,  and  ZIP  Code)  10  SOURCE  OF  FUNDING  NUMBERS 

RollinP  Air  PnrrA  Raoe.  PROGRAM  I  PROJECT  TASK  (WORK  UNIT 

Oiling  Air  rorce  Base  element  no  no  no  accession  n< 

Washington  DC  20332  ,  r 

_ _ _ _ _  61102F  2307  ° _ 1 

1 1  TITLE  (Include  Security  Classification) 

A  Numerical  Simulation  of  the  Flow  Field  and  Heat  Transfer  in  a 
Rectangular  Passage  with  a  Turbulence  Promoter 

12  PERSONAL  AUTHOR(S)  “  .  ““ 

Bryan  R.  Becker  and  Richard  B.  Rivir 

13a  TYPE  OF  REPORT  3b  TIME  COVERED  ''  Jl4.  DATE  OF  REPORT  (tear,  Month,  Day)  |15  PAGE  COUNT 

FINAL  FROM_Jun_87_  to  Sep  87  Bine  19 PA  c>? 

16  SUPPLEMENTARY  NOTATION 

1 7. _ ’  _ CQSATI  CODES _  18  SUBJ^CT^Tf  RMS  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

_ FIELDI  /  GROUP  SU8-GKOUP  1 

_ 1  02 _ Heat  Transfer,  Turbulence  Augmentation  -  )  ■ 

19  ABSTft(^T  ( Continue  on  reverse  if  necessary  and  identify  by  block  number)  . —  ••••-■•  . . 

A  numerical  study  of  the  flow  field  and  heat  transfer  <r  a  turbine  blade  intern?! 
cooling  passage  is  described.  Attention  is  focused  cr.  a  short  straight  section  of  the 
rectangular  passage  winch  includes  a  single  turbulator  protruding  from  its  floor.  The 
tvr  dimensional,  transient,  Reynolds  averaged  Navier  Stokes,  contiruiJy  and  energy 
equations  are  iterated  to  a  steady  state  solution  using  the  MacConnack  explicit  predictor- 
corrector  algorithm,  lurbulence  closure  is  achieved  through  the  use  of  the  Baldwin- 
Lomax  forr  of  the  Cebe*.  1 -Smith  algebraic  tw^  layer  eddy  vi’scositv  model. 

Plots  of  skin  friction,  local  heat  transfer'  rate,  streamlines,  velocity  profiles  and 
temperature  profiles  are  given.  It  war  found  that  the  widely  used  Reynolds  Analogy 
greatly  underpredicts  the  heat  transfer  rale  as  qiven  by  a  direct  calculation  usirq 
Fourier's  lav/.  _ _ _ _ _ 


17,  1 

COSATI  CODES 

FIELDI 

/  GROUP  SU8-GKOUP 

11  1 

/  92 

1  04 

20  distribution /av'"  't'/ of  eg.,  «act 

Kj  UNCLASSIFiED'UNLiM!  TED  □  SAME  AS  RPT 
22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 

RICHARD  B.  RIVIR 

DD  Form  1473,  JUN  86 


?  1  ARSTRACT  TfCURiT'-  CLASSIFICATION 

□  one,  users  UNCLASSIFIED _ _ 

22b  TELEPHONE  (Include  Area  Code,  |  ,','c  OFF-Ct  SYvyiOL 

513/255-6768  !  AFWAL/POTC 

F'eviow  editions  aie  obsolete  _ SI  CURIVY  CLASSlFIC  A* ION  CV  T  his  ‘  A^f 

UNCLASSIFIED 


rl/»/  *  »*  <,*  4_*  »*  »/ 


TABLE  OF  CONTENTS 


SECTION  PACE 

>  I  Introduction .  1 

i  II  The  Numerical  Model .  7 

A.  Computational  Domain . 7 

,  B.  Governing  Equations .  8 

i  C.  Boundary  Conditions .  11 

: 

» 

D.  Initial  Conditions .  14 

.  E.  Numerical  Algorithm .  15 

i 

|  III  Results .  16 

■  A.  Convergence .  16 

■  B.  Velocity  and  Temperature  Profiles .  18 

p  C.  Streamlines,  Skin  Friction,  and  Heat  Transfer 

|  Rates .  22 

1  IV  Conclusions .  28 

i 

(References .  30 

Nomenclature .  32 


FIGURE 


LIST  OF  ILLUSTRATIONS 


K 


I 


I 


1  The  100  by  60,  (X,Y),  computational  mesh  used  in 

this  numeiical  study . 

2  Horizontal  velocity  profile  at  K  =  1 . 

3  Temperature  profile  at  K  =  1 . 

4  Vertical  velocity  profile  at  K  =  1 . 

5  Horizontal  velocity  profile  at  K  =  6 . 

6  Temperature  profile  at  K  =  6 . 

7  Vertical  velocity  prolile  at  K  =  6 . 

8  Horizontal  velocity  profile  at  K  =  11 . 

9  Temperature  profile  at  K  =  11 . 

10  Vertical  velocity  profile  at  K  =  11 . 

11  Horizontal  velocity  profile  at  K  =  18 . 

12  Temperature  profile  at  K  =  18 . 

13  Vertical  velocity  profile  at  K  =  18 . 

14  Horizontal  velocity  profile  at  K  -  2d . 

15  Temperature  profile  at  K  =  78 . 

16  Vertical  velocity  profile  at  K  -  28 . 

17  Horizontal  velocity  prolile  at  K  =  S3 . 

18  Temperature  profile  at  K  =  33 . 

19  Vertical  velocity  prolile  at  K  =  33 . 

20  Horizontal  velocity  profile  at  K  -  38 . 

31  Temperature  profile  at  K  =  38 . 

22  Vertical  velocity  profile  at  K  -  38 . 

23  Horizontal  velocity  profile  at  k  -  47 . 

24  Temperature  profile  at  r  -  <7 . 

25  Vertical  velocity  profile  at  k  -  4 . 

28  Horizontal  volcci."  plot  le  at  k  =  w . 


passages.  Such  prediction  is  necessary  to  achieve  a  thermal  design  of 
these  components  which  incorporates  efficient  cooling  and  thereby 
alleviates  the  effects  of  the  increased  heat  load. 


There  have  been  numerous  experimental  and  analytical  studies  of  the 
fluid  flow  and  heat  transfer  within  straight  rectangular  ducts  which 
simulate  the  flow  within  the  internal  cooling  passages  inside  a  turbine 
blade.  Several  investigators  have  studied  the  phenomena  of  secondary 
flows  in  ducts  with  noncircular  cross  section.  These  secondary  flows, 
transverse  to  the  main  flow  direction,  are  composed  of  eight  counter¬ 
rotating  vortices  which  convect  momentum,  vorticity  and  total  energy 
into  the  corners  of  the  duct. 

in  an  early  study,  Gessner  and  Jones  [6]  utilized  hot  wire 
anemometry  to  study  secondary  flow  effects  in  fully  developed  turbulent 
air  flows  within  rectangular  channels.  They  studied  flows  in  channels 
with  smooth  interna]  surfaces  and  reported  that  the  secondary  flow 
velocity,  when  normalized  bv  the  axial  mean  flow  velocity,  decreased 
with  increasing  Reynolds  number.  They  also  deduced  that  secondary 
flows,  composed  of  eight  counterrotating  streamwise  vortices,  are  caused 
by  an  imbalance  between  Reynolds  stresses  and  static  pressure  gradients. 
However,  they  presented  no  heat  transfer  results. 

Launder  and  Ying  [ LO |  expanded  upon  the  work  of  Cessner  and  Jones 
[6]  by  again  using  hot  wire  anemometry  to  study  secondary  air  flows  in  a 
duct  of  square  cross  section.  However,  C bey  considered  ducts  whose 
surfaces  were  roughened  with  transverse  mounted  square  turbulators  and 
concluded  that  the  secondary  flow  velocity,  when  normalized  by  the 
average  friction  velocity,  became  independent  of  Reynolds  numbei  and 
surface  roughness.  Again.,  no  heat  transfer  results  were  reported. 
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Later,  Gessner  [5]  reported  a  detailed  analysis  of  the  mechanisms 
which  initiate  secondary  flow  in  developing  turbulent  air  flow  along  a 
corner.  He  concluded  that  transverse  gradients  of  the  Reynolds  shear 
stress  components  are  responsible  for  the  generation  of  secondary  flow 
within  turbulent  flows  along  a  corner.  Also,  he  reiterated  that 
secondary  flows  act  as  a  significant  convective  transport  agent  of  the 
momentum,  vorticity  and  total  energy  of  the  main  axial  motion. 

Melling  and  Whitelaw  [12]  were  the  first  investigators  to  apply 
laser-Doppler  anemometry  to  the  study  of  secondary  flows  within 
developing  turbulent  flow  inside  rectangular  ducts.  They  used  water  as 
the  working  fluid  and  their  results  confirm  the  qualitative  results 
published  by  earlier  investigators  but  they  indicated  quantitative 
differences.  Their  data  are  valuable  for  the  verification  of 
mathematical  turbulence  models.  Also,  their  work  showed  that  the 
developing  flow  is  much  more  sensitive  to  inlet  conditions  than  is  a 
nominally  fully  developed  flow.  They  did  not  study  heat  transfer 
effects . 

In  summary,  the  focus  of  the  above  mentioned  experimental 
investigations  was  upon  the  cause  and  structure  of  secondary  flows 
within  turbulent  flows  inside  rectangular  ducts.  None  of  these 
inves t 1 ga t ions  studied  the  heat  transfer  phenomena  occurring  within  these 
duct  flows. 

Experimental  studies  more  closely  related  to  the  thermal  design  of 
turbine  blade  cooling  passages  began  to  appear  in  the  1980' s.  For 
example,  Han  [7]  studied  the  effects  of  square  lurbulator  geometry  upon 
the  heat  transfer  coefficient  and  friction  factor  in  fully  developed 
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turbulent  air  flow  through  a  square  duct  with  two  opposite  ribbed  walls. 


He  used  a  blower  to  force  room  air  through  a  test  section  consisting  of 
four  heated  aluminum  plates  of  which  the  two  vertical  plates  were 
ribbed.  The  heaters  were  independently  controlled  so  as  to  provide  a 
controllable  constant  heat  flux  throughout  the  test  section.  The  test 
section  was  instrumented  with  thermocouples  and  pressure  taps  while  the 
inlet  air  flowrate  and  temperature  as  well  as  the  exit  air  temperature 
were  monitored.  The  friction  factor  was  determined  from  the  pressure 
drop  data  and  the  Stanton  number  was  determined  from  the  air  and  wall 
temperatures  with  the  constant  wall  heat  flux  as  a  knowrn  boundary 
condition.  In  this  investigation  Han  studied  the  effects  of  rib  height 
to  equivalent  duct  diameter  (H/D)  and  rib  height  to  rib  pitch  or  spacing 
along  the  duct  (H/P) .  He  reports  that  both  the  friction  factor  and  the 
Stanton  number  increase  with  increasing  H/D  and  likewise  for  increasing 
H/P.  He  also  presents  an  empirical  prediction  method  for  determining 
the  friction  factor  and  the  Stanton  number  for  flows  with  this  geometry. 

In  an  additional  study  using  the  same  experimental  techniques,  Han 
et  al  [9]  investigated  the  combined  effects  of  rib  angle  of  attack  (u) 
and  rib  height  to  pitch  (H/P).  They  report  that  the  maximum  heat 
transfer  and  friction  factor  occurred  at  a  of  60-75  degrees.  They  also 
report  that  as  a  decreased  from  90  degrees  to  approximately  45-30 
degrees,  • ne  heat  conductance,  at  a  constant  pressure  drop,  increased  to 
its  maximum  value  while  the  pressure  diop  required  for  a  constant  heat 
transfer  rate  decreased  to  its  minimum  value.  Thus,  it  appears  that  an 
angle  of  attack  in  the  range  of  45-30  degrees  produces  the  optimal 
thermal  performance  lor  a  given  ax  in.  pressure  drop.  They  also  report 
that  as  tlie  H/P  ratio  decreases  the  same  trends  are  evident  but  are 


relatively  reduced. 


Finally,  Han  et  al  1.8]  measured  the  combined  effects  of  rib  angle 
of  attack  (a)  and  duct  aspect  ratio.  Whereas  previous  work  by  Han  had 
been  concerned  with  ducts  of  square  cross  section,  this  study  involved 
rectangular  ducts  whose  ribbed  walls  were  2  to  4  times  wider  than  their 
smooth  walls.  They  found  that,  for  this  type  of  duct,  the  heat  transfer 
and  friction  factor  varied  only  slightly  as  a  varied  from  45  degrees  to 
90  degrees,  with  the  maximums  occurring  at  a  equal  to  90  degrees.  They 
also  found  that  the  heat  transfer  and  friction  factor  increased  as  the 
width  of  the  ribbed  wall  increased  relative  to  the  width  of  the  smooth 
wall . 

Several  numerical  studies  of  the  flows  in  rectangular  channels  with 
turbulators  have  also  been  reported.  For  example,  Durst  and  Rastogi 
[3]  performed  both  experimental  and  numerical  studies  of  the  turbulent 
flow  over  a  transverse  mounted  square  turbulator  in  a  two  dimensional 
channel.  Their  experimental  work  involved  the  use  ot  a  one  dimensional 
laser-Doppler  anemometer  to  measure  turbulence  characteristics  in  a 
plane  two  dimensional  water  channel .  They  used  the  TtACH  numerical 
algorithm  with  a  K-e  turbulence  model  to  simulate  the  two  dimensional 
flow  in  the  recirculation  region  near  the  turbulator.  However,  the 
velocity  ami  turbulence  profiles  near  the  wall  as  well  as  the  wall  shear 
stress  anu  hence  the  heat  transfer  at  the  wall  were  all  determined  by  a 
"law  of  the  wall"  rather  than  calculated  as  part  of  a  detailed  numerical 
analysis . 

Fujita  et  al  [4]  numerically  analyzed  fully  developed  turbulent 
flow  in  a  square  duct  witli  two  opposite  roughened  walls.  They  used  a 
stream  func tion-vorticity  formulation  to  solve  for  the  flow  field  in 
pianes  orthogonal  to  the  main  f Low  direction  anu  then  marched  this 


solution  in  the  streamwise  direction.  Their  results  show  the  presence 


of  only  four  streamwise  vortices  rather  than  the  eight  reported  by  the 
experimentalists  mentioned  above.  Also,  their  calculations  were 
performed  on  a  very  coarse  11  by  11  uniform  mesh  over  a  quadrant  of  the 
cross  section  of  a  rectangular  duct  and  hence  do  rot  capture  much  of  the 
detail  of  the  heat  transfer  phenomena  near  the  solid  wall. 

Linton  and  Shang  [llj  investigated  jet  impingement  cooling  of 
turbine  blades  in  a  rotating  frame  of  reference.  Tney  numerically 
solved  the  laminar  three-dimensional,  compressible  Navier-Stokes 
equations  to  yield  a  description  of  the  flow  field  for  both  rotating  and 
stationary  cases.  However,  their  calculations  omitted  both  the  effects 
of  turbulence  and  surface  roughness. 

In  conclusion,  the  authors  feel  that  there  remains  a  need  for  the 
detailed  numerical  modeling  of  the  heat  transfer  and  fluid  flow 
processes  very  near  the  fluid-solid  interface.  This  paper  reports  on  an 
initial,  two  dimensional  modeling  effort  which  begins  to  address  this 
need.  In  terms  of  long  range  goals,  this  detailed  numerical  model  will 
be  expanded  to  include  the  effects  of  rotation,  multiple  turbulence 
promoters  and  three  dimensional  secondary  flows.  In  this  way,  as  more 
detail  and  completeness  is  added  to  the  model,  the  local  heat  transfer 
rate  along  the  cooling  passage  will  be  more  accurately  predicted, 
thereby  enabling  the  optimization  <>f  cooling  configuration  design  prior 
to  full  scale  prototype  testing  and  at  a  greatly  reduced  cost.  This 
optimal  heat  transfer  design  will  result  in  lower  metal  temperatures  and 
reduced  thermal  gradients,  thus  enabling  the  turbine  Slades  to  withstand 
the  increased  operating  conditions  necessary  to  improve  engine  efficiency 
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SECTION  II 


THE  NUMERICAL  MODEL 


A.  Computational  Domain 

In  this  initial  numerical  study  of  the  fluid  llow  and  heat  transfer 
in  a  turbine  blade  internal  cooling  passage,  attention  is  focused  on  a 
two  dimensional  flow  over  a  single  turbulator  of  approximately  square 
cross  section  which  protrudes  upward  from  the  floor  of  the  passage, 
transverse  to  the  main  flow  direction.  The  X-coordinate  is  oriented 
along  the  length  of  the  passage,  coincident  with  the  bulk  flow  direction 
while  the  Y-coordinate  is  oriented  in  the  direction  of  the  turbulator 
height,  normal  to  tfie  passage  floor.  This  two  dimensional 
approximation,  which  greatly  reduces  the  computational  complexity  and 
cost,  is  considered  to  be  sufficient  for  this  initial  study.  The 
accuracy  of  this  two  dimensional  assumption  increases  as  the  passage 
aspect  ratio  (height/transverse  width)  decreases  and  is  a  very 
i  asonabie  numerical  approximation  to  the  experimental  study  reported  by 
Han  et  al  [bj  for  a  passage  aspect  of  0.25. 

An  outline  of  the  computational  mesh  or  domain  is  shown  in  Figure 
I.  The  main  flow  direction  is  from  left  to  right  with  the  left  and 
right  boundaries  being  the  inflow  and  outflow  boundaries,  respectively. 
The  lower  boundary  is  a  solid  wall  boundary  while  the  top  boundary  is  a 
symmetry  boundary  which  simulates  the  centerline  of  the  flow  passage. 

The  turbulence  promoter  or  rib  is  modelled  as  an  obstruction  protruding 
upward  from  the  bottom  boundary  witb  a  flat  top  whose  length  in  the 
streamwisc  direction  is  equal  to  its  height.  The  upstream  and  down¬ 
stream  feces  of  the  rib  re  modelled  as  sin  [ fx  •  0/(2  •  rib  height)]. 
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so  that  the  total  width  of  the  rib  is  three  rib  heights,  measured  from 
the  point  of  departure  from  horizontal  to  the  point  of  return  to 
horizontal . 

The  computational  domain  extends  20.0  rib  heights  in  the  streamwise 
X-direction  end  4.2  rib  heights  in  the  vertical  Y-direction  with  the 
upstream  face  of  the  rib  beginning  its  departure  from  horizontal  at 
approximately  4.0  rib  heights  from  the  inflow  boundary.  There  are  100 
grid  lines  in  the  X-direction,  clustered  around  the  turbulence  promoter 
with  the  distance  between  grid  lines  growing  geometrically  from  an 
initial  spacing  of  0.1  rib  heights.  There  are  60  grid  lines  in  the 
vertical  Y-direction  with  44  within  the  first  rib  height  which  is  well 
within  the  inrbulent  boundary  layer.  The  mesh  spacing  grows  in  the 
Y-direction  as  a  geometric  progression  from  an  initial  spacing  of  0.01 
rib  heights.  From  this  description,  it  can  be  seen  that  the 
computational  mesh  has  very  high  resolution  and  is  capable  of  revealing 
considerable  detail  of  the  flow  field  ard  heat  transfer  near  the 
turbulence  promoter. 


B.  Governing  Equations 

The  governing  equations  for  this  simulation  consist  of  the  two 
dimension.!,  Reynolds  averaged  form  of  the  following  turbulent, 
compress!'  ii..  unsteady  equations  -  continuity,  viscous  Navier  Stokes, 
and  energ; : 


->  v  •  (on)  =  0 
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The  dependent  variables  in  these  equations  are  the  time  averaged  values 


of  density,  x -momentum,  y-momentum  and  total  energy:  p,  pu,  pv,  and 

pe  =  p [cvT  +  i  (u2  +  v2) J . 

The  Reynolds  averaging  of  the  turbulent  equations  produces 
components  of  Reynolds  stress  which  appear  in  the  stress  tensor,  t,  and 
correlations  of  velocity  and  temperature  fluctuations  which  appear  in 
the  heat  flux  vector,  q.  The  required  turbulence  closure  is  achieved 
through  the  use  of  the  Baldwin-Lomax  form  of  the  Cebeci-Smith  algebraic 
two  layer  eddy  viscosity  model  [1],  The  Reynolds  stress  components  are 
then  modelled  as  the  product  of  the  eddy  viscosity,  t,  times  the 
velocity  gradient  of  the  mean  flow  while  the  correlations  of  velocity 
and  temperature  fluctuations  are  represented  with  a  turbulent  Prandtl 
number  (Pr(  0.9): 

3u.  3u.  „  3U.  //N 
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,P  ,  e  \  3T  (5) 
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Sutherland's  law  is  used  to  calculate  the  molecular  viscosity  of 
air,  u,  at;  a  function  of  temperature.  The  molecular  Prandtl  number,  Pr, 
is  assumed  to  have  a  constant  ’'alue  of  0.73  while  the  constant  volume 
specific  heat  of  air,  c  .  is  taken  to  be  4c90.0  ft2/s2-°R:  (>  =  c  /c 

V  P  V 

1.4).  Th.-  equation  of  state  is  assumed  to  be  the  ideal  gas  law  with  the 
gas  constant  for  air  equal  to  17ib.U  i  t  ‘7::‘’-°R. 

As  si  "Wn  in  Figure  1,  the  (X,Y)  cooidinate  system  is  a  body  fitted 
coordinate  system  which  is  deformed  to  account  for  the  rib  on  the  lower 
boundary.  This  deformed  coord -nut e  system  results  in  computational 
cells  which  are  not  rectai-gr  lai  bur  'allies  distorted  quadrilaterals.  In 


i 
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order  to  maintain  second  order  calculational  accuracy  in  the  spatial 


domain,  a  coordinate  transformation  is  made  to  a  uniform  rectangular 
(£,n)  grid.  by  means  of  the  chain  rule,  the  governing  equations  in  the 
transformed  space  can  be  written  in  conservation  law  form  as  follows: 
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In  Equation  (6),  £ 
derivatives  of  the 
to  the  coordinates 
follows: 


3?  .  ,  3(5  .  n  (6) 

nx  ^  +  S  H  \  =  ° 

,  n  »  £  »  and  q  are  the  first  order  partial 
x  x  y  y 

transformed  independent  variables  (£,n)  with  respect 
(X,Y)  ,  and  the  vectors  fj ,  f  and  2  are  defined  as 


(7) 


(8) 


(9) 


This  set  of  equations  along  with  the  appropriate  boundary  and  initial 
conditions  constitutes  the  mathematical  model  of  the  cooling  air  flow 
through  a  rectangular  passage  with  a  turbulence  promoter. 


SK 


C.  Boundary  Conditions 

Specification  of  appropriate  and  consistent  boundary  conditions  is 
essential  to  achieving  a  realistic  and  stable  numerical  solution  for  any 
fluid  flow  or  heat  transler  problem.  In  the  current  simulation  there 
are  four  boundaries  which  must  be  considered:  bottom,  top,  inflow,  and 
outflow.  Since  there  are  four  unknowns,  then  at  each  boundary,  four 
conditions  must  be  specified. 

The  bottom  boundary,  J  -  1,  is  a  solid,  isothermal,  no  slip 
boundary  maintained  at  T  ^  =  580. 0°R,  which  implies  the  following 
conditions : 

(Pu)J=1  =  (Pv)J=1  =  0  (10) 

(pe)  .  =  p  ..c  T  „  (11) 

J=1  wall  v  wall 

The  assumption  that  the  normal  derivative  of  the  pressure  vanishes  at 
the  wall,  coupled  with  the  specification  of  the  wall  temperature, 
implies  the  following  condition  on  the  density  at  the  wall: 

‘  “wall  '  "2) 

wall 

The  top  boundary,  J  ~  JMAX,  which  is  iocated  at  the  centerline  of 
the  rectangular  passage,  is  a  symmetry  boundary  within  the  free  stream. 
Symmetry  implies  that  the  vertical  velocity  vanishes: 

(PV).’MAX  =  0 

Since  the  top  boundary  is  within  the  free  stream,  the  horizontal 
velocity,  pressure,  temperature,  and  density  ere  all  equal  to  the  free 
stream  vciues  which  are  specified  as  follows: 


u  =  262.21  ft/s 

CO 

P  =2116.8  lb, /ft2 

°°  t 

T  =  540. 0°R 


P  =  0.002284  lb  *s2/f r 4 
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These  specifications  imply  the  following  boundary  conditions: 

PJMAX  =  P~  (18) 

(PU)JMAX  =  P~U«  (19) 

(pe)JMAX  "  P~  fCvT^+  *  (U~)2]  (20) 

At  the  upstream  inflow  boundary,  I  =  i,  it  is  assumed  that  the  flow 

experiences  reversible,  adiabatic  acceleration  from  stagnation  conditions 
in  a  settling  tank  upstream  of  the  inflow  boundary.  To  determine  the 
corresponding  stagnation  temperature  and  pressure,  it  is  assumed  that 
the  air  in  the  settling  tank  contains  the  total  energy  of  the  free 
stream.  That  is,  the  stagnation  temperature  and  pressure  are  determined 
by  applying  the  first  law  of  thermodynamics  and  Bernoulli's  equation  to 
a  reversible  adiabatic  total  deceleration  from  the  free  stream  velocity, 
pressure,  and  temperature  specified  above  in  Equations  (14)  through 
(16).  The  resulting  stagnation  temperature  and  pressure  are  thus 
calculated  to  have  the  following  values: 

T  =  545 . 72°R  (21) 

o 

Pq  =  2195.33  lbf / ft2  (22) 

These  upstream  settling  tank  conditions  are  then  used  to  determine  the 
inflow  boundary  conditions.  First,  it  is  assumed  that  the  flow  has  been 
conditioned  prior  to  entering  the  rectangular  passage  so  that  it  enters 
as  parallel  flow  with  no  vertical  component  whicli  imposes  the  following 
boundary  condition  on  vertical  momentum: 

(pv)I=1  =  0  (23) 

Second,  it  is  assumed  that  the  normal  gradient  of  the  temperature 
vanishes  at  the  inflow  which  implies  that  die  temperature  is 
extrapolated  trom  the  interior  of  the  flow  field: 


PW 


T  =  i 
1=1  1=2 


Thus,  the  horizontal  velocity  component  can  be  calculated  by  applying 
the  first  law  to  the  reversible  adiabatic  acceleration  of  the  air  from 
the  settling  tank: 

“i-i  -  'VVTi»1/2  <“> 

The  isentropic  pressure  relation  is  then  used  to  determine  the  pressure 


at  the  inflow: 


pi-i '  WV 


Y/Y-l 


Finally,  the  boundary  condition  on  the  density  at  the  inflow  is  given  by 
using  the  ideal  gas  law: 

pi=i  “  pi/,(rairti)  (27) 


(pu)l=1  =  pjUj  (28) 

(pe)I=1  “  PitcvTi  +  i  (29) 

At  the  downstream  outflow  boundary,  I  =  IMAX,  the  pressure  is 
specified  to  be  the  free  stream  value  while  the  normal  gradients  of  the 
velocities  and  total  energy  are  set  to  zero  which  implies  that  the 
values  of  these  variables  are  extrapolated  from  the  interior  of  the  flow 


field : 


P  =  P 
IMAX 


UIMAX  UIMAX-1 


VIMAX  V1MAX-1 


CIMA/  elMAX-l 


The  conditions  on  the  velocities  and  the  total  energy  at  the  outflow 
imply  a  spc rif ication  of  the  temperature  at  the  outflow  which  in  turn 
implies  a  boundary  condition  on  the  density  at  the  outflow  via  the  ideal 
gas  law: 


PIMAX  =  PIMAX  1  (RAIR*TIMAX) 


^PU^IMAX  =  PIMAXUIMAX 
<‘PVjIMAX  =  PIMAXVIMAX 
(‘pe)IMAX  =  °lMAXeIMAX 


D.  Initial  Conditions 

In  the  current  simulation,  no  attempt  was  made  to  produce  a  time 
accurate  solution  of  the  initial  transient  of  the  developing  flow. 
Therefore,  the  initial  conditions  were  used  only  as  a  starting  point 
from  which  the  algorithm  could  begin  calculation  of  a  steady  state 
solution. 

To  this  end,  at  the  beginning  of  the  calculation,  it  was  assumed 
that  a  turbulent  boundary  layer,  two  rib  heights  in  thickness,  6, 
existed  throughout  the  rectangular  channel.  Within  this  boundary  layer, 
the  horizontal  velocity  was  given  by  the  1/7  power  law: 

u(y)  =  ^^(y/S)1^7  (38) 

Above  the  boundary  layer,  the  horizontal  velocity  was  set  equal  to  the 
free  stream  velocity.  Equation  (14). 

Throughout  the  flow  field,  the  vertical  component  of  velocity  was 
set  to  zero;  while  the  pressure,  temperature  and  density  were  set  equal 
to  the  free  stream  values,  Equations  (15)  -  (17).  The  initial 
condition::  on  the  momentum  and  total  energy  were  then  calculated  using 
the  initial  values  of  these  primitive  variables.  Finally,  the  density 
and  energy  at  the  solid  boundary  were  modified  to  account  for  the  higher 
temperature  of  the  wall. 
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E.  Numerical  Algorithm 

The  current  mathematical  model  of  the  cooling  air  flow  through  a 
rectangular  passage  with  a  turbulence  promoter  is  composed  of  the 
computational  domain,  the  set  of  governing  equations,  the  boundary 
conditions,  and  the  initial  conditions.  The  numerical  algorithm  used  to 
calculate  the  solution  to  this  model  is  the  MacCormack  explicit 
predictor-corrector  technique  as  implemented  by  Shang  and  Hankey  [14]. 

Although  this  computer  code  can  be  utilized  to  achieve  a  time 
accurate  transient  solution,  it  also  contains  an  option  which  allows 
accelerated  convergence  to  a  steady  state  solution  through  the  use  of 
local  time  stepping.  This  local  time  stepping  option  was  invoked  in  the 


SECTION  III 


RESULTS 


A.  Convergence 

In  Shang's  implementation  of  the  MacCormack  algorithm,  various 
options  exist  to  accelerate  convergence  and/or  maintain  numerical 
stability.  The  algorithm  can  be  run  in  the  local  time  step  mode  to 
accelerate  convergence  or  the  time  accurate  mode  to  maintain  stability. 
The  algorithm  contains  a  numerical  smoothing,  which  when  invoked 
promotes  stability,  but  can  also  degrade  the  solution.  Finally,  the 
turbulent  viscosity  can  be  omitted  from  the  iteration  which  promotes 
stability  but  also  produces  a  less  realistic  solution.  Therefore,  a 
strategy  must  be  developed  for  the  use  of  these  various  options  so  as  to 
achieve  a  stable  but  realistic  solution  while  minimizing  the  number  of 
iterations  required.  For  example,  the  iteration  could  be  initiated  with 
all  options  turned  to  full  stability.  Then,  after  a  number  of 
iterations,  options  could  be  switched  to  maximize  the  rate  of 
convergence.  Finally,  a  number  of  iterations  should  be  run  with  the 
options  switched  to  produce  the  most  realistic  solution.  The  strategy 
used  in  tie  current  simulation  is  shown  in  Table  1. 

In  Snang's  implementation,  convergence  is  measured  by  the  L-2  norm 
of  the  change  in  the  normalized  global  solution  vector  between 
successive  iterations,  ||6L'n+j|J.  The  global  solution  vector,  U,  is  the 
concatenation  of  all  of  the  nodal  solution  vectors  as  given  in  Equation 
(7).  Each  component  of  the  change  in  the  global  solution  vector  is 
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Iteration 

Time 

Smoothing 

Turbulent 

1 

Count 

Step 

Viscosity 

0-  5,000 

5,001-10,000 

10,001-45,000 

45,001-99,999 


Time 

Accurate 

Local 
Time  Step 

Local 
Time  Step 
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squared,  normalized  by  the  free  stream  value  of  the  respective  variable 


and  summed: 
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Typically,  in  the  iteration  to  a  converged  steady  state  solution,  the 
value  of  this  norm  will  drop  3  orders  of  magnitude  from  the  initial 
conditions  and  then  oscillate  about  this  reduced  value,  indicating  that 
convergence  has  been  achieved  [13].  As  shown  in  Table  2  for  the  current 
simulation,  the  value  of  this  norm  decreased  from  7.137  at  the  outset  to 
0.002662  after  50,000  iterations  and  oscillated  near  that  value  for  the 
remaining  49,999  iterations.  This  behavior  is  typical  of  convergence  to 
the  steady  state  solution.  The  current  simulation  required  92  minutes 
of  execution  time  on  a  Cray-1  computer. 


B.  Velocity  and  Temperature  Profiles 

As  described  in  Section  II. A.,  the  computational  grid  is  20  rib 
heights  in  length  with  the  top  front  corner  of  the  rib  located  5  rib 
heights  downstream  from  the  inflow  boundary,  and  the  top  back,  corner  of 
the  rib  located  6  rib  heights  downstream  from  the  inflow  boundary. 

There  are  100  grid  lines  in  the  streamwise  X-direction,  clustered  around 
the  rib.  These  grid  lines  are  numbered  as  the  K  index  starting  at  the 
inflow  boundary.  Hence,  K  =  1  corresponds  to  the  inflow  boundary. 

Profiles  of  u/u  ,  v/v  ,  and  T/T  were  plotted  against  Y/YMAX  at 
the  16  streamwise  locations  listed  in  Table  3.  In  these  profile 
plots,  given  in  Figures  2  through  49,  a  rib  height  corresponds  to 
Y/YMAX  =  0.237.  Hence,  in  the  profile  plots  for  K  values  of  28,  33,  38, 


the  curve  does  not  reach  Y/YMAX  =  0  because  the  bottom  of  the 


computational  mesh  coincides  with  the  top  of  the  rib  or  turbulence 
promoter  which  has  a  Y/YMAX  value  greater  than  zero. 

The  profiles  at  K  =  1  represent  the  inflow  boundary  conditions 
where  the  flow  has  not  yet  properly  developed.  This  explains  the 
uniformly  zero  vertical  velocity  and  the  fluctuating  horizontal 
velocity. 

At  K  *  6,  the  profiles  show  chat  the  momentum  and  thermal  boundary 
layers  are  beginning  to  develop.  Since  this  flow  is  predominately  in 
the  horizontal  direction,  the  vertical  velocity  is  very  small  compared 

> 

» 

|  to  the  horizontal  velocity.  Therefore,  the  small  variations,  which  are 

* 

>  due  to  round  off  in  the  calculation  of  the  vertical  velocity,  are  nearly 

i 

* 

j  the  same  magnitude  as  the  actual  vertical  velocity,  thus  producing  the 

exaggerated  fluctuation  exhibited  In  the  vertical  velocity  profile. 

'  At  K  *  11,  and  18,  the  profiles  show  the  continued  development  of 

"  the  momentum  and  thermal  boundary  layers  and  a  vertical  acceleration  as 

j  the  flow  begins  to  move  up  and  over  the  turbulence  promoter.  Also  at  K 

►  =  18,  both  the  vertical  and  horizontal  velocities  are  negative  near  the 

solid  wall  which  indicates  the  presence  of  a  recirculation  zone  which 
will  be  discussed  later. 

The  locations,  K  =  28.  33,  and  38,  are  at  the  front  top  corner, 
midline  and  back  top  corner  of  the  rib,  respectively.  The  vertical 
velocity  is  positive,  indicating  that  the  floe,-  is  moving  up  over  the 
rib.  The  horizontal  velocity  is  greatly  accelerated  near  the  solid 
wall,  exhibiting  a  uniform  proiiie  at  K  =  38,  typical  of  a  fully 
developed  turbulent  I  low.  Both  the  momentum  and  therma  l  boundary  layers 
have  been  destroyed  cue  to  the  intense  turbulent  roving  which  is 


occurring.  The  temperature  near  the  wall  is  much  lower  than  at  previous 
streamwise  locations  and  exhibits  some  numerical  overshoot.  This 
overshoot  as  well  as  the  fluctuations  in  both  temperature  and  horizontal 
velocity  profiles  are  numerical  artifacts  produced  by  over  restrictive 
boundary  conditions  at  the  top  boundary  which  will  be  corrected  in 
future  calculations. 

The  streamwise  locations,  K  =  47  through  K  =  83,  sta’t  at  the  back 
bottom  corner  of  the  rib  and  extend  to  the  reattachment  point  which 
marks  the  end  of  the  large  recirculation  bubble  downwind  of  the  rib. 

The  temperature  profiles  show  that  the  turbulent  mixing  behind  the  rib 
has  greatly  diffused  the  energy  from  the  hot  solid  wall  to  form  a  warm 
air  region  which  extends  almost  2  rib  heights  vertically  from  the  solid 
wall.  The  horizontal  velocity  profile  is  typical  of  a  recirculation 
zone  with  a  negative  upstream  flow  near  the  wall  changing  to  a  positive 
downstream  flow  higher  in  the  flow  field.  The  vertical  velocity 
profiles  also  show  evidence  of  flow  recirculation. 

The  profiles  at  K  =  85,  89,  and  100,  uemonstrate  that  both  the 
momentum  and  thermal  boundary  layers  are  reestablished  downstream  of  the 
recirculation  bubble.  However,  comparison  of  these  profiles  to  those  at 
K  =  b,  11,  and  18,  upstream  of  the  turbulence  pronofer,  show  that  the 
downstream  boundary  layers  are  much  thicker  and  well  mixed,  illustrating 
the  effects  of  the  intense  turbulent  mixing  caused  by  the  flow  over  the 


turbulence  promoter. 


Streamlines,  Skin  Friction,  and  Heat  Transfer  Rates 


Careful  analysis  of  the  printed  output  at  each  of  the  100 
streamwise  locations  reveals  the  interesting  details  of  the  flow  pattern 
shown  in  Figure  50.  There  are  five  recirculation  zones  or  bubbles 
produced  by  the  flow  over  the  turbulence  promoter.  There  is  a  clockwise 
rotating  zone,  1.0  rib  height  in  length,  which  begins  1.0  rib  height 
upstream  of  the  rib  and  is  approximately  0.1  rib  height  tall.  The 
second  clockwise  rotating  bubble,  which  is  about  a  half  rib  height  in 
length  and  very  thin,  begins  just  downstream  from  the  top  front  corner 
of  the  rib.  Behind  the  rib  there  are  two  counterc lockwi se  vortices 
which  are  separated  by  a  clockwise  vortex.  The  first  clockwise  vortex, 
which  is  1.0  rib  height  in  length,  begins  just  behind  the  back  top 
corner  and  is  0.54  rib  height  in  thickness.  Adjacent  to  this  clockwise 
vortex  is  a  counterclockwise  vortex  of  approximately  the  same  size. 
Finally,  there  is  the  very  large  recirculation  bubble  which  begins  just 
downstream  of  the  back  top  corner  of  the  rib  and  attains  a  maximum 
thickness  of  0.71  rib  height  just  behind  the  rib.  This  large  bubble 
extends  to  the  reattachment  point  which  is  located  approximately  7.0  rib 
heights  from  the  back  top  corner  or  6.0  rib  heights  from  the  back  bottom 
corner . 

The  location  of  tiie  reattachment  point  agrees  very  well  with  the 
wind  tunnel  investigations  of  flow  over  a  surface  mounted  cube  conducted 
by  Castro  and  Robins  [ 2  j .  They  report  that  the  wake  produced  by  a 
uniform  turbulent  stream  completely  decayed  within  b.O  cube  heights 


downstream. 


Skin  friction  and  heat  transfer  rates  were  calculated  at  the  16 


streamwise  locations  listed  in  Table  3.  The  values  of  local  skin 
friction,  shown  in  Table  4  and  plotted  in  Figure  51,  were  calculated 
from  the  following  equation: 

,  Twall  UV\=0  (40) 


The  normal  derivative  of  the  horizontal  velocity  was  approximated  as  a 
central  difference  between  the  J  =  1  node  (solid  wall)  and  the  J  =  2 
node  (Y  <  0.3281  E-4  feet).  The  viscosity  was  determined  by  averaging 
the  values  given  by  Sutherland's  Law  evaluated  at  Tj_^  and  Tj_„. 
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Table  4.  Local  Skin  Friction  Coefficient 


K 

(Streamwi.se 

Index) 


Skin  Friction  Coefficient 
(c|  x  104) 
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The  local  heat  transfer  rates,  given  in  Table  5  and  plotted  in 
Figure  52,  are  expressed  interms  of  the  nondimensional  Stanton  number: 


CpP»U“(Twall~T“) 

These  heat  transfer  rates  were  calculated  by  two  different  methods.  In 
the  first  method,  the  Reynolds  Analogy  was  used,  in  which  the  heat 
transfer  rate  is  related  to  the  skin  friction  as  follows: 

Stanton  Number  =  ic^  ^2) 

This  method  is  widely  used  in  turbomachinery  calculations  but  these 
calculations  greatly  underpredict  the  heat  transfer  rates  measured 
experimentally . 

In  the  second  method,  Fourier's  Conduction  Law  was  used  to  evaluate 
the  heat  transfer  rate  at  the  wall,  Qwa^]_»  which  appears  in  Equation 


0  -  (43) 
Qwall  M3yV0 

In  this  calculation,  the  thermal  conductivity  was  treated  as  a  constant 
and  the  normal  derivative  of  temperature  was  evaluated  in  a  manner 
similar  to  that  used  for  the  normal  derivative  of  the  horizontal 
velocity  discussed  above. 

In  the  convection  heat  transfer  process,  the  energy  is  first 
transferred  by  conduction  from  the  solio  wall  to  the  fluid  particles 
adjacent  to  the  wall.  These  fluid  particles  are  then  transported  away 
from  the  solid  wall  by  the  bulk  mixing  of  the  fluid  which  results  in  the 
diffusion  of  energy  from  the  solid  wall.  The  basis  for  the  Reynolds 
Analogy  is  the  assumption  that  the  same  mechanism  is  responsible  for  the 
exchange  of  both  heat  and  momentum.  This  assumption  disregards  the 
details  of  the  energy  transfer  process  at  the  solid  wall-fluid 


interface.  As  shown  in  Table  5  and  Figure  52,  the  direct  calculation  of 
the  heat  transfer  rate  via  Fourier's  Law  predicts  heat  transfer  rates 
which  can  be  as  great  as  30  tines  that  given  by  the  Reynolds  Analogy. 
This  great  discrepancy  can  be  explained  by  the  fact  that  the  Reynolds 
Analogy  is  derived  for  laminar  or  fully  turbulent  flows  over  a  flat 
plate  with  laminar  and  turbulent  Prandtl  number,  Pr  =  Pr  =  1.0.  Thus, 
the  Reynolds  Analogy  is  not  applicable  to  flows  with  recirculation  zones 
and  flow  reversals  which  cause  the  normal  gradient  of  velocity  to  differ 
greatly  from  the  normal  gradient  of  temperature. 


SECTION  IV 


CONCLUSIONS 


The  flow  field  and  heat  tran&Ier  along  a  turbine  blade  internal 
cooling  passage  was  numerically  simulated  using  a  two  dimensional 
geometry.  Attention  was  focused  on  the  flow  along  a  short  straight 
section  of  a  rectangular  passage  with  a  single  transverse  mounted 
turbulator  protruding  from  its  floor.  The  two  dimensional,  transient, 
Reynolds  averaged  Navier  Stokes,  continuity  and  energy  equations  were 
iterated  to  a  steady  state  solution  using  the  second  order  accurate 
MacCormack  explicit  predictor-corrector  algorithm.  Turbulence  closure 
was  achieved  through  the  use  of  the  Baldwin-Lomax  form  of  the 
Cebeci-Smith  algebraic  two  layer  eddy  viscosity  model. 

In  contrast  to  the  work  done  by  earlier  investigators,  this  study 
was  performed  using  a  computational  mesh  which  had  very  high  resolution. 
There  were  100  grid  lines  in  the  horizontal  (axial)  direction  clustered 
around  the  turbulator  and  60  grid  lines  in  the  vertical  direction  with 
44  within  the  first  rib  height  which  was  well  within  the  turbulent 
boundary  layer.  This  high  resolution  mesh  revealed  considerable  detail 
of  the  flow  field  and  heat  transfer  near  the  turbulator. 

The  calculation  revealed  a  complex  flow  structure  near  the 
turbulator  with  five  stationary  or  standing  vortices.  This  flow 
structure,  which  is  more  complex  than  is  ordinarily  expected  for  such  a 
simple  passage/turbulator  geometry,  has  a  pronounced  effect  upon  the 
heat  transfer  as  shown  m  Table  5  and  Figure  52.  This  significant 
insight  indicates  that  lurtiier  study  is  needed  into  che  effects  of 
turbulator  geometry  and  spacing  (or  pitch)  a.:  a  means  cf  heat  transfer 


enhancement . 


Furthermore,  this  calculation  revealed  a  considerable  difference 
between  the  heat  transfer  rate  calculated  by  using  the  Reynolds  Analogy 
and  that  given  by  a  direct  calculation  using  Fourier's  law.  This 
discrepancy  results  from  the  complex  flow  structure  near  the  turbulator 
which  causes  the  normal  gradient  of  velocity  to  differ  greatly  from  the 
normal  gradient  of  temperature.  Similar  flow  structures  exist 
throughout  many  turbomachinery  components,  and  therefore,  corresponding 
heat  transfer  effects  may  occur  in  those  components  which  would  also  be 
inadequately  predicted  by  the  Reynolds  Analogy.  This  significant 
insight  may  lead  to  remedies  for  the  underprediction  of  heat  transfer 
rates  in  many  turbomachinery  calculations. 
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NOMENCLATURE 


Cj  =  skin  friction 

0^  =  constant  pressure  specific  heat  of  air 

c^  =  constant  volume  specific  heat  of  air 

D  =  equivalent  passage  diameter 

e  =  total  energy 

F  =  flux  vector 

G  ■  flux  vector 

H  *  rib  or  turbulator  height 

I  =  horizontal  index 

IMAX  *  maximum  horizontal  index 

T  =  vertical  index 

JMAX  =  maximum  vertical  index 

P  =  pitch,  spacing  between  ribs  or  turbulators 
=  flow  field  pressure 


Pr 

=  molecular 

Prandtl 

number 

Prt 

=  turbulent 

Prandcl 

number 

Q.q 

-  heat  flux 

St 

=  Stanton  number 

T  =  air  temperature 

t  =  t ime 

-► 

U  =  solution  vector 


u  =  horizontal  velocity 

X, x  =  horizontal,  axial  coordinate 

Y, y  =  vertical  coordinate 

YMAX  =  maximum  vertical  coordinate 


-  *  •  *  m.w  -  *  , 
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NOMENCLATURE  (Continued) 


Greek 

a  =  rib  or  turbulator  angle  of  attack 

Y  =  ratio  of  specific  heats 

6  =  turbulent  boundary  layer  thickness 

6.,  =  Kronecker  delta 

ij 

e  =  turbulent  eddy  viscosity 

n  =  transformed  coordinate 

U  “  molecular  viscosity  of  air 

£  =  transformed  coordinate 

p  =  air  density 

i  =  stress  tensor 

Subscript 

o  »  total  values 

co  ~  free  stream  value 

wall  =  value  at  solid  boundary 


>«r 


xvmA 


o 
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Reynolds  Analogy 


